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WIENER CHAOS SOLUTIONS OF LINEAR STOCHASTIC 
EVOLUTION EQUATIONS 

BY S. V. LOTOTSKY 1 AND B. L. ROZOVSKII 2 

University of Southern California 

A new method is described for constructing a generalized solution 
of a stochastic evolution equation. Existence, uniqueness, regularity 
and a probabilistic representation of this Wiener Chaos solution are 
established for a large class of equations. As an application of the gen- 
eral theory, new results are obtained for several types of the passive 
scalar equation. 

1. Introduction. Consider a stochastic evolution equation 

(1.1) du{t) = {Au{t) + f{t)) dt + (Mu(t) + g(t)) dW(t), 

where A and Ai are differential operators, and W is a cylindrical Brownian 
motion on a stochastic basis F= ($l,J-, {J r t}t>o^)- Let Aq and .Mo be the 
leading (highest-order) terms of the operators A and J\4, respectively, and 

Z t (u)= f{Mu{s)+g(s))dW(s). 
Jo 

Traditionally, (1.1) is studied under the following assumptions: 

(i) The operator Aq-\MqM* is elliptic. 

(ii) The noise term Zt(v) is sufficiently regular. More specifically, for a 
suitable function space X, E J T ||-2t(u)||x dt < oo for all v in a dense subspace 
of X. 

Under these assumptions, there exists a unique Ito (strong) solution or a 
martingale (weak) solution u of (1.1) so that u € L2^l X (0,T);X) for T > 
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(see, e.g., [5, 19, 33, 35]). In the future, we will refer to such solutions as 
traditional or square integrable solutions. 

There are important examples demonstrating that the assumptions (i)- 
(ii) are necessary for the existence of a square integrable solution of (1.1). 

In particular, it was shown in [30] that the stochastic advection-diffusion 
equation 

d 

— u(t, x) = Au(t, x) + u(t, x)W(t, x), x E R d , 

where W(t,x) is the space-time white noise, has no square integrable solu- 
tions in any Sobolev space X = H s (M. d ), s 6l, if d > 2. In this case assump- 
tion (i) holds but assumption (ii) does not. 
Simple calculations show (see, e.g., [35]) that 

d 2 d 
du(t, x) = a 2 ——^u(t, x) dt + a—u(t, x) dw(t), 
ox* ox 

where w(t) is a one-dimensional Brownian motion and u(0,x) is square in- 
tegrable, has no square integrable solutions unless a 2 — \o~ 2 > 0, that is, 
unless (i) holds. Obviously, (i) implies that the order of the operator M 
must be no larger than a half of the order of A. Moreover, this assumption 
is, in some way, counterintuitive. Indeed, in the deterministic theory, only 
the highest-order operator (.Ao, in our case) has to be elliptic. 

The objective of the current paper is to study stochastic differential equa- 
tions of the type (1.1) without the restrictive assumptions (i) and (ii). The 
basic idea can be described as follows. If (1.1) does have a sufficiently reg- 
ular solution, this solution can be projected on an orthonormal basis in 
some Hilbert space, resulting in a system of equations for the corresponding 
Fourier coefficients. We now turn this argument around, and define the so- 
lution of (1.1) as a formal Fourier series with the coefficients computed by 
solving the corresponding system. It often happens that this system has a 
solution under more general conditions than the original equation. 

This approach could be traced to the classical separation of variables 
ideas in PDE's. For example, the Navier-Stokes equation is often defined as 
a system of coupled ODE's for the modes of its formal Fourier expansion 
with respect to the spatial variables (see, e.g., [7] and [26]). Similarly, in 
our case, the nonrandom spatiotemporal variables (x,t) are being separated 
from the "random variable" (Brownian motion). 

More specifically, the traditional solution of (1.1) is an Tjf -measurable 
square integrable random variable taking values in a Hilbert space X. A 
classical result by Cameron and Martin [4] provides an orthonormal basis 
3 = {£,a-,ct G J} in the Wiener Chaos space L2(W;X), where J is the set 
of multi-indices a = {af} of finite length \a\ = J2i°k=i a i- Accordingly, a 
Wiener Chaos solution of (1.1) is defined as a formal Fourier series with 
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respect to the Cameron-Martin basis H. By construction, this solution is 
strong in the probabilistic sense, that is, uniquely determined by the coeffi- 
cients, free terms, initial condition and the Wiener process. The coefficients 
in the Fourier series are computed by solving the corresponding propagator, 
a lower-triangular system of deterministic parabolic equations, uniquely de- 
termined by (1.1) (in earlier works, e.g., [28], the propagator is referred to 
as S-system). 

Of course, unless both (i) and (ii) hold, one could not expect the resulting 
Fourier series to converge in L2(W;X). However, we demonstrate that un- 
der quite general assumptions, the Fourier series converges in a "minimal" 
weighted Wiener Chaos space L2,q(W;X). The construction of this space is 
quite simple (cf. [30]). Given a sequence of positive numbers Q = {qi,q2, ■ ■ ■}, 
we define L2,q(W;X) as the collection of sequences u = {u a ,a E J} with 
u a G X so that 

II n 2 let 1 1 ii 2 „ 

IMIl 2 ,q(W;X) : = 2^ 1 IKHx<°0> 

where 

q a = Y[qf% a£j. 

i,k 

We remark that the space L2(W; X) is a very special example of the sequence 
spaces studied in [15, 16, 17]. 

Of course, the Wiener Chaos solution is a weak solution designed for deal- 
ing with equations that do not have square integrable solutions. However, in 
some problems, the Wiener Chaos solution serves as a convenient first step 
in the investigation of square integrable solutions. For example, this is the 
case in the proof of the existence of square integrable solutions of degenerate 
parabolic SPDE's (Corollary 4.2 and Theorem 6.3). 

Constructions based on various forms of the Wiener Chaos decomposition 
are popular in the study of stochastic differential equations, both ordinary 
and with partial derivatives. For stochastic ordinary differential equations, 
[20] used multiple Wiener integral expansion to study Ito's diffusions with 
nonsmooth coefficients. More recently, LeJan and Raimond [23] used a sim- 
ilar approach in the construction of stochastic flows. Various versions of the 
Wiener Chaos appear in a number of papers on nonlinear filtering and re- 
lated topics; see, for example, [3, 24, 28, 31, 37]. The book by Holden et al. 
[12] presents a systematic approach to the stochastic differential equations 
based on the white noise theory; see also [11, 34] and the references therein. 

The propagator was first introduced by Mikulevicius and Rozovskii [27], 
and further studied in [24] , as a numerical tool for solving the Zakai filtering 
equation. In [30], the propagator was used to construct a generalized solu- 
tion of the reaction-diffusion equation driven by the space-time white noise 
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in several space dimensions. A similar system can be derived for certain 
nonlinear equations, such as the stochastic Navier— Stokes equation [29]. 
The main results of the paper are: 

(1) Existence, uniqueness, regularity and the Krylov-Veretennikov for- 
mula for the Wiener Chaos solution of (1.1) (Theorems 3.4, 4.1 and Corol- 
lary 4.2). 

(2) A Feynman-Kac formula for Wiener Chaos solutions in L2,q(W; )) 
(Theorem 5.1). 

(3) Existence, uniqueness and regularity properties for the transport equa- 
tion with: 

(a) a space-time white noise- type velocity field (Theorem 6.1); 

(b) an incompressible Kraichnan turbulent velocity field, the nonvis- 
cous case (Theorem 6.3). 

The paper is organized as follows. Section 2 introduces the Cameron- 
Martin basis and the weighted Wiener Chaos spaces. Section 3 presents the 
general definition of the Wiener Chaos solution and establishes the con- 
nection with the traditional and the white noise solutions. In the three 
main Sections, 4, 5 and 6, three types of stochastic equations are stud- 
ied under three different sets of assumptions; these assumptions are listed 
at the beginning of the corresponding section. Section 4 presents the basic 
existence/uniqueness/regularity result for an abstract stochastic evolution 
equation under assumptions (Al)-(A3). Section 5 establishes probabilistic 
and multiple Wiener integral representations of Wiener Chaos solutions for 
stochastic partial differential equations in M. d with nonrandom coefficients 
under assumptions (B0)-(B5). Finally, Section 6 illustrates the results of the 
previous two sections for several versions of the turbulent transport equation 
under assumptions (S1)-(S3). 

The following notation will be in force throughout the paper. A is the 
Laplace operator, Di = d/dxi, i = 1, . . . ,d, and summation over the repeated 
indices is assumed. The space of continuous functions is denoted by C, and 
H](R d ), 7 G R, is the Sobolev space 

if- l/(y)| 2 (l + |y| 2 ) 7 dy < 00 >, where / is the Fourier transform of /. 
I m d J 

2. Weighted Wiener Chaos spaces. In this section we review the con- 
struction of the Cameron-Martin basis and define the spaces of generalized 
random elements. 

For a fixed T > 0, let F = (fi, {^ r t}o<t<T,I I> ) be a stochastic basis with 
the usual assumptions and let W = (wk = ttffc(i), k > 1, < t < T) be a col- 
lection of independent standard Brownian motions on F. Denote by Tf 
the (T-algebra generated by the random variables {wk(s),k > l,s< t}, and 
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by L2(W), the Hilbert space of Jj! -measurable square integrable random 
variables. 

Our first step is to construct the Cameron-Martin basis, a special or- 
thonormal basis in the space L2(W). 

Let m = {rrik,k > 1} be an orthonormal basis in L2((0,T)) so that each 
function rrik = mfc(t) is bounded for t G [0,T]. Given such a basis m, define 
independent standard Gaussian random variables 

f T 

Cik = / rrii{s)dw k {s). 
Jo 

Consider the set of multi-indices 

J= la = (a$,i,k> l),af e {0, 1,2, .. < oo L 

I i,k ) 

The set J is countable, and, for every a G JT, only finitely many of a- 3 are 
not equal to zero. For a G JT, we write 

l«l=I]af, a!=II a ^ ! ' 
and define the collection S = {£ a ,a G ,7} of random variables so that 
(2-1) £« = 4rII 

where 

if„(t) = (-l)"e* 2 /^ e -* 2 /2 

is nth Hermite polynomial. 

Theorem 2.1. TTie collection H = {^ a ,a G ^7} is an orthonormal basis 
in L 2 (W). 

Proof. This is a version of the classical result by Cameron and Martin 
[4]. □ 

Using the Cameron-Martin basis E in L2(W), we now define the space of 
generalized random elements. 

Let X be a Banach space and Q = {qi,q2, . ■ ■}, a sequence of positive 
numbers. Define 

(2-2) <? a = IKS 
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Definition 2.2. (i) The Q-weighted Wiener Chaos space L 2 ,q(W;X) 
is the collection of sequences u = {u a ,a G J} with u a G X so that 

II n 2 2a 1 1 n 2 ^ 

IMIl 2 ,q(W;X) : = ^ IK|| X <°°- 

(ii) The Q~-weighted Wiener Chaos space L 2j q- (W; X) is the collection of 
sequences u = {u a ,a G ^7} with u a G X so that 

II l|2 _ -2ai| ||2 

IpIIl 2Q _(w ; x) : ~ ^ II u qIIx <00 - 

The space X will be omitted from the notation if X = R. The sequence 
Q will be omitted from the notation if Q = 1, that is, if q k = 1 for all k. 

Given a u = {ii a ,a G J} G L2,q(W;X), we call each u a a generalized 
Fourier coefficient of it and identify u with a formal Fourier series 

The members of the set L2 i q(W;X) are called X-wa/uec/ generalized random 
elements. Similarly, the members of the set L2,q(W;Li((0, T);X)) are called 
'X.-valued generalized random processes. 

For u G L 2 ,q(W; X) and v G L 2>Q - ( W ) we define 

(2.3) = XI UaVa > 

the series in (2.3) converges in the norm of X by the Cauchy-Schwartz 
inequality. 

3. Linear stochastic evolution equations and the propagator. In this sec- 
tion we define the Wiener Chaos solution and establish its connection with 
the traditional and white noise solutions. To motivate the definition, we start 
by reviewing the main results about the traditional solution. 

Let (V, H, V') be a normal triple of Hilbert spaces so that V C H C V 
with both embeddings continuous; for the complete definition of the normal 
triple see Section 3.1 in [35]. Denote by (v',v), v' G V, v G V, the duality 
between V and V' relative to the inner product in H. 

For t G [0, T], consider families of linear operators A = A(t) and = 
Mk(t) so that, for each t, the operators A(t) : V — ► V , Mk{t) : V — > H are 
bounded. Consider the following equation: 

(3.1) = u o + j {Au{s) + f(s)) ds + J (M k u(s) + g k (s)) dw k {s), 

< t < T. 
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Recall that summation convention over the pairs of repeated indices is in 
force. 

We proceed with a review of the traditional approach. Assume that, for 
a31v€V,t€ [0,T], 

(3.2) £ll>l*(t>>l&<°o, 

k>\ 

and the nonrandom input data uq , / and g k satisfy 

rT r T 



(3.3) \\u f H +[ \\f(t)\\ 2 v>dt + J2 [ 

Jo k>1 Jo 



\9k(t)\\ H dt < °°. 



Definition 3.1. An J^-adapted process u £ L 2 (W;L 2 {(0,T);V)) is 
called a square integrable solution of (3.1) if, for every v G V, there exists 
a measurable subset Q' of with P(fl') = 1, so that, for all < t < T, the 
equality 

(u(t),v) H = {u ,v) H + / (Au(s) + f(s),v) ds 

(3.4) t J ° 

+ Y1 / (M k u(s) + gk{s),v) H dw k (s) 
k>i Jo 

holds on W. Similarly, u G Li2((0,T);V) is a solution of the deterministic 
equation 

u(t) = uq + / Au(s)ds+ / f(s)ds 
Jo Jo 

if, for every v G V and £ G [0,T], equality (3.4) holds with M. k = g k = 0. 

Existence and uniqueness of the traditional solution of (3.1) are estab- 
lished under an additional assumption about the operators A and M. k . 

Definition 3.2. Equation (3.1) is called strongly parabolic if there ex- 
ists a positive number e and a real number Co so that, for all v G V and 

te[0,T], 

(3.5) 2(X(t)u,i;) + ]T ||Mfc(£H|^ + e|M|^> < C \\vf H . 

k>l 

Equation (3.1) is called weakly parabolic (or degenerate parabolic) if condi- 
tion (3.5) holds with e = 0. 



8 S. V. LOTOTSKY AND B. L. ROZOVSKII 

Theorem 3.3. If (3.3) holds and (3.1) is strongly parabolic, then there 
exists a unique square integrable solution of (3.1). The solution process u 
belongs to 

£ 2 (W; L 2 ((0, T)-V)) n L 2 (W; C((0, T),H)) 

and satisfies 



\ 2 v dt 



(3.6) 



E( sup \\u(t)f H + f \\u(t)\\ 

\0<t<T JO 

<C(C ,e,T)(\\u \\ 2 H + f T \\f( t )\\ 2 v/ dt+J2 f T \\9k 



%dt 



Proof. This follows, for example, from Theorem 3.1.4 in [35]. □ 

When (3.1) is weakly parabolic, then the solvability result is somewhat 
different; see Section 3.2 in [35] for details. 

As an element of the Hilbert space L 2 (W; L2((0,T); V)), the traditional 
solution of (3.1) admits a representation u(t) = J2 a ej u aif)ia hi the Cameron- 
Martin basis H. 

Theorem 3.4. An T% -adapted process u is a square integrable solution 
of (3.1) if and only if u(t) = Yl<a£j u a{t)£,a so that the Fourier coefficients 
u a satisfy 

(3.7) ( I* \\u a {t)\\ 2 v dt+ sup IK^lH) <oo, 

a( zj\J0 0<t<T / 

and solve the propagator 

u a (t) = u I(\a\ = 0) + / (Au a (s) + f(s)I(\a\ =0))ds 



(3.8) 

+ J ^\[<4( M kU a -^ k ){s) +g k (s)I(\a\ = l))mi(s)ds, 



/0 i, k 

where a~(i,k) is the multi-index with components 



(a (i,k)) 



! j max(a t - —1,0), if i = j and k = I, 



a 1 - 



Before presenting the proof, we define the Wiener Chaos solution of (3.1). 
The definition is motivated by Theorem 3.4. 
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Definition 3.5. A V- valued generalized random process u is called a 
Wiener Chaos solution of (3.1) if the generalized Fourier coefficients u a ,a G 
J , of u are a solution of the propagator (3.8). 

To prove Theorem 3.4 and to derive an alternative characterization of the 
Wiener Chaos solution, we need a few additional constructions. 
Denote by Ti. the set 

H= (J Loo((0,T);M n ). 
n>l 

If h G 7i, then there exists an JV > 1 so that 

h= (hi,...,h N ), 
with each fo^ G Loo((0,T)). We define 



f(t,/i) = exp^(j[ 



fc fc (s)du> fc (a)-§ j |/i fe (s)| 2 ds) ), 



(3.9) £(h)=£(T,h), 

hik = / h k (t)rrii(t) dt, mi G m, 
Jo 

and 

The following properties of the process £(t,h) are verified by direct calcula- 
tion: 

h a 



(3.10) £{h) = J2 



'CC\ 



(3.11) £{t,h) = l+ [ £(s,h)h k (s)dw k (s). 

Jo 

One consequence of (3.10) is that £(h) G L2,q(W) for every weight sequence 
Q, with ll«?(^)ll| 2 Q (w) = ex P(EfcLi9fcl|/ifc|ll 2 ((o,T)))- Another consequence of 
(3.10) is an alternative representation of £ a : 



(3.i2) ^7SW m 



h=0 



Remark 3.6. By Lemma 4.3.2 in [32], the family {£(h),h G H) is dense 
in L2(W) and therefore in every iv2Q(W). 
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Proof of Theorem 3.4. (i) Assume that u = u(t) is a traditional 
solution of (3.1). Equation (3.12) implies 

ua(t) = E(«(t)&) = -^=^E(u(t)£(h)) 



h=0 

Using the .F^-measur ability of u(t) and the martingale property (3.11) of 
£(t,h), we derive 

E(u(t)£(h))=E(u(t)E(£(h)\Fl v ))=E(u(t)£(t,h)), 

and (3.8) follows after applying the Ito formula to the product u(t)£(t,h) 
and differentiating the resulting equation with respect to h. By Theorem 
3.3.3(iii) in [22], this differentiation with respect to h is justified. 
The same arguments show that the time evolution of 

(3.13) ut)-=m a \^) 

is described by 

(3.14) £ a (t)=I{\a\=0)+ fY.^t^-^k){s)m(s)dw k {s). 

Jo i,k 

(ii) Conversely assume that condition (3.7) holds. Then the process u(t) = 
J2aej v>a(t)€a satisfies 

u £ L 2 (W; L 2 ((0, T); V)) n L 2 (W; C((0, T);H)), 

and the function u h =E(u£(h)) belongs to L 2 ((0,T); V) flC((0,T);iJ). By 
Parsevall's identity and relation (3.10), 

x , h a 

uh= y. —^u a . 

Using (3.8), we conclude that, for every ip € V and h € T~C, the function Uh 
satisfies 



(u h (t),(p) H = (u ,(p) H + J (Au h (s),(f) H ds+ I (f(s),<p)ds 
h c 



h r i — 



a£j ' i,k 



+ (9k(s),<p) H I(\a\ = l))ds. 



To simplify the above equality note that if I(t) = Jq (M k u(s), ip)n dwk(s), 
then the Fj^-measurability of I(t) and relation (3.14) imply 



(3.15) E(I(t)£ a )=E(I(t)Ut)) = J o Y^^ m i( s )(M k u a - m (s),<p) H ds 
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Similarly, 
El 



ia ( gk(s)dw k (s)] =V / Ja%mi(s)(g k (s),<p) H I(]a\ = l)ds. 



Therefore, 

h c 
a! 

a£j i,k 



h f I — 

77 / V a i m i( s )(( M kU a ~( l ,k)(s),^>) H + (9k(s),<f) H I(\a\ =l))ds 







As a result, 



E(£ <p) H ) = E(£(h)(u , <p) H ) +E[S(h) [ (Au(s), <p) ds 



o 



(3.16) +E(s(h) j\f(s),<p)ds\ 



+ K\S(h) ((Mku(s),<p) H + (g k (s),<p) H )dw k (s) 

Equality (3.16) and Remark 3.6 imply that, for each t and each ip, (3.4) 
holds with probability 1. Due to continuity of u, a single probability-1 set 
can be chosen for all t € [0, T] . 
Theorem 3.4 is proved. □ 

Theorem 3.4 implies that, if it exists, the traditional solution of (3.1) 
coincides with the Wiener Chaos solution. 

We now give an alternative characterization of the Wiener Chaos solution. 



Theorem 3.7. A V -valued generalized random process u is a Wiener 
Chaos solution of (3.1) if and only if, for every h&TC, the function u^{t) = 
{(u(t) , £ (h))) is a solution of the equation 

(3.17) u h {t)=u + [ (Au h (s) + f(s) + h k (s)M k u h {s) + h k (s)g k (s))ds. 
Jo 

Proof. Assume that u is a Wiener Chaos solution of (3.1), that is, 
u £ L2,q(W; V) for some Q and each u a is a solution of (3.8). By definition 
(2.3) and relation (3.10), 



(3.18) 
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Then (3.17) follows from (3.8). Indeed, 

Yj ~n Y V a i M k U a-(i.k)m = X] Yj I -I- ,^ M k u a -(i,k)mihik 

( h a \ 

= X] zJ —J= M kU a mih ik = h k M k u h . 



Lk \ a J(X 



Computations for the other terms are similar. 

Conversely, assume that u € L2,q(W; ^((0, T); V)) and Uh € ^((0, T);V) 
is a solution of (3.17). By relation (3.18), 



1 <9l Q l 



\dh a 



h=0 



Then term-by-term differentiation in (3.17) implies (3.8). Theorem 3.7 is 
proved. □ 

To conclude this section, we briefly discuss the relation between the 
Wiener Chaos and white noise solutions. While the white noise solution 
for (3.1) has not been defined in general, analysis of the particular cases 
[11, 12, 34] shows that the white noise solution exists provided (3.17) has a 
solution that is an analytic function of h. The white noise solution of (3.1) 
is constructed as an element of the space 

(<S)_ P (X) = < u a : X r a/ll M «llx < 00 f° r some £ > >, 



where 



i ■ k 

(S)-o(M. d ) is known as the space of Hida distributions [11]. Even though 
(<S)_ i0 (X) is not, in general, related to L2,q(W;X), in many examples one 
can take either q k = q or q k = Cq k for some q < 1. For such weights, 

£ 2 ,q(W;X)c(cS)_ (X) 

with strict inclusion, so that the Wiener Chaos approach provides better reg- 
ularity results. In addition, we remark that in contrast to the Wiener Chaos 
solutions, the white noise solutions are weak solutions in the probabilistic 
sense. 
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4. Linear evolution equations in weighted Wiener Chaos spaces. In this 
section we present the main existence, uniqueness and regularity result for 
the Wiener Chaos solution of (3.1). We make the following assumptions: 

(Al) There exist positive numbers C\ and 5 so that, for every v £ V and 

t£[0,T], 

(4.1) (A(t)v,v)+8\\v\$<C x \\vf H . 

(A2) There exists a real number C*2 and a sequence of positive numbers 
Q = {ik, k > 1} so that, for every v £ H and t £ [0, T], 

(4.2) 2(A(t)v,v) + ]T q 2 k \\M k (t)v\\ 2 H < C 2 \\v\\ 2 H . 

k>l 

(A3) The initial condition uq is nonrandom and belongs to H; the process 
/ = f(t) is deterministic and J \\f(t)\\y, dt < 00; each g\. = g k (t) is a 
deterministic process and 



E / Ak(t)\\ 2 H dt< 00. 



Denote by (Pj, s ,0 < s <t <T) the semigroup generated by the operator 
A; Pt := Ptfl- In other words, the solution of the equation 

v (t) = v + [ Av{r)dT+ I f(r)dT, < s < t < T, 

J S J S 

v £H,fe L 2 ((0,T);V), is written as 

v(t) = P t>s v + j* P t , T f(r)dT. 

This solution exists by assumption (Al) and Theorem 3.3. Denote by ur \ 
the solution of (3.8) when \a\ =0. 

Theorem 4.1. Under assumptions (A1)-(A3), (3.1) has a unique Wiener 
Chaos solution. The solution u = u(t) has the following properties: 

(1) For every 7 G (0,1), 

(4.3) u £ L 2i7Q (W; L 2 ((0, T); V)) n L 2nQ (W; C((0, T);H)). 

(2) For every < t < T, u(t) £ L2 i q(W; H) and the following relations 
hold: 

\a\=n 

= E /"*" 

k\ , . . . , k n ^ 1 



Jo 
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(4-4) ••• r Pt,s n M kn ---P S2tSl 







x (M kl U(o) + ?fciflfci(si)) rfwfci(si) • • • dw kn (s n ), 

n > 1, 



w/iere .Mfc = q k M k ; 

\l 2 ,q(W;H) 



\<t)\\ 2 l 



(45) <3e<**(||tio||2r + C/ jf* 11/0011^ * + £ jTflgllflfcWIIir^, 

where the number C2 is from (4.2) and the positive number Ct depends only 
on 5 and C\ from (4.1). 

Proof. The arguments build on the techniques from [25]. 

(1) By Theorem 3.3, there exists a unique traditional solution of equation 

v(t)= u + f {Av{s) + f(s))ds 







(4.6) 

+ £ / (M k v(s) + g k (s))jq k dw k (s). 
k>i Jo 

By Theorems 3.4 and 3.7, v a = r )\ a \q a u a , and (4.3) follows. 

(2) Direct computations show that, for \a\ = 1 with a* = 1, the corre- 
sponding solution u a := ur^ of (3.8) is 

(4.7) U(i k ^(t)=J P t , s (MkU( )(s) + g k (s))mi(s)ds. 

Since &i k \ = £i k = Jq rrii(s) dw k (s), we conclude that 

^Uik(t)C( ik ) = I Pt, s {M k U(Q)(s) + g k (s)) dw k (s) 



i>l 

or 



q a u a (t)£ a = Y^ I Pt,s{MkU(Q){s) +q k g k (s))dw k (s). 
\a\=l k>l J ° 

Continuing by induction on \a\ and using the relation between the Hermite 
polynomials and the iterated ltd integrals ([13], Theorem 3.1), we derive 
(4.4). 
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An immediate consequence of (4.4) is the following energy equality: 

E <? 2Q IK(*)lll 

\a\=n 

(4-8) = E 



fci,...,fe„>i 



o jo 



r s 2 2 

J \\Pt,s n M kn ■ ■ ■ P S2 , Sl {M kl u {0) + q kl g kl )\\ H ds n , 



where ds n = ds\- • ■ ds n . 

To derive (4.5), consider first the homogeneous equation with / = g k = 
and define F n (t) = '^2\ a i =n q 2a \\u a (t)\\'fj, n > 0. Note that u^(t) = P t u . 

By assumption (4.2), 

(4.9) j t F (t) < C 2 F (t) - E \\M k P t u \\ 2 H . 

k>l 



For n > 1, the energy equality (4.8) implies 

^n(i)= E / / •••/ 

at , . , jo jo Jo 



ft rs n pS2 

(4.10) + E / / •••/ 2(APt, Sn M kn ---M kl P Sl u , 
k u ...,k n >i Jo Jo Jo 

Pt,s n M kn ■ ■■M kl P Sl u ) ds n . 



By assumption (4.2) 

rt rs n rS2 



ft rs n rS2 

E //•••/ 2(AP t , Sn M kn ---M kl P Sl u , 
kl ,...,k n >i Jo Jo Jo 

(4.11) Pt, Sn M kn • • • M kl P sl u ) ds n 



ft fS„ fS2 

k ^ JO JO JO 



fcl,...,fc n >l 
As a result, for n > 1 

at 



. /■< fSn-l fS2 _ 

(4.12)+ E / / •••/ \\M kn P t , Sn _ 1 M kn _ 1 ---M kl P sl u \\ 2 H ds n - 

,. i -.Jo JO JO 



fcl,...,fc n >l 
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rt rs n rS2 



E / H--- n\\Mk n+1 Pt,s n M kn ---M kl P Sl u f H ds n . 
^ Jo Jo Jo 



fcl,...,fe„ + l>l 
Consequently, 

N N 

dt 



( 4 - 13 ) i E E i 2a h*{t)f H < c 2 E E <i 2a h«m 2 H, 

n=0\a\=n n=0 \a\=n 



so that, by the Gronwall inequality, 

N 

(4.14) £ J2 q 2a \\n a (t)f H <e c ^\\u f H 

n=0 |a|=n 

or 

(4-15) \Ht)\\l 2Mw . H) <e c * t \\ Uo \\ 2 H . 

The remaining two cases, namely, no = gk = and ito = / = 0, are analyzed 
in the same way, and then (4.5) follows by the triangle inequality. Theorem 
4.1 is proved. □ 

Corollary 4.2. Let (Hj^h^^ik^k be deterministic measurable func- 
tions of (t,x) so that 

\aij(t,x)\ + \bi(t,x)\ + \c(t,x)\ + \a ik (t,x)\ + \u k (t,x)\ < K, 

i, j = 1,. .. ,d,k > l,x G R d ,0 <t < T; 

(%•(*, x) - \a ik {t,x)a jk (t,x))y i y j > 0, 

x,yeR d ,0<t<T; and 

E \vk(t,x)\ 2 < C v < oo, 
fc>i 

x€lR rf ,0<t<r. Consider the equation 

du = (Di(aijDju) + biDiU + cu + f) dt 

(4.16) 

+ (a ik DiU + z; fc u + g k ) dw k . 

Assume that the input data satisfy u G L 2 (M d ) ; / G L 2 ((0,T); i7^ 1 (R d )), 
Sfc>i ll5felli 2 ((o T)xR d ) ^ °°' anc ^ there exists an e > so i/taf 

a ij (t,x)y l y j >e\y\ 2 , x, y G M d , < t < T. 

Then there exists a unique Wiener Chaos solution u = u(t,x) of (4.16). The 
solution has the following regularity: 

(4.17) u(V) £L 2 (W;L 2 (R d )), 0<t<T, 
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and 



(4.18) 



Ellu 




LaCflVT);^ 0*")) 



E ll5 f fellL 2 ((0,T)xIR d ) 
k>l 



where the positive number C* depends only on C U ,K,T and e. 



Proof. Direct computation shows that the operators 



A = D i (a ij D j ) + b i D i + c, 



M k = o- ik Di + v k 



satisfy assumptions (Al) and (A2) in the normal triple {H^iW*) , L2(M d ), 
if 2 ~ 1 (M d )) with q k = 1. Then relations (4.17) and (4.18) follow from (4.5). 

□ 

Note that in contrast to the statement of this corollary, all previous results 
concerning (4.16) required additional regularity of the coefficients and input 
data; see, for example, [35], Section 4.2. 

Theorem 4.1 is a bona fide extension of Theorem 3.3. Indeed, if condi- 
tion (3.5) holds so that (3.1) is strongly parabolic, then, taking Q = 1, we 
recover the statement of Theorem 3.3. Further analysis of condition (3.5) 
indicates that, for a strongly parabolic equation, one can find an admissible 
weight sequence Q so that q k = q > 1 . Since for Q > 1 we have a strict in- 
clusion L2,q(W;X) C L2(W;X), Theorem 4.1 represents an improvement of 
Theorem 3.3 for strongly parabolic equations. 

For weakly parabolic equations, similarly to the proof of Corollary 4.2, 
one can take Q = 1, and then Theorem 4.1 represents an extension of the 
existing results ([35], Sections 3.2 and 4.2). 

For nonparabolic equations, the results of Theorem 4.1 are completely 
new. 

Remark 4.3. The formal Fourier series u = J2asj u a£,a for the Wiener 
Chaos solution is a generalization of the representation formula for solutions 
of SODE's, derived by ([20], Theorem 4) using iterated Ito integrals; see also 
[23], Theorem 3.2. Indeed, 



and equality (4.4) connects the inner sum on the right with the iterated 
integrals. See also Example 5.2 below. 
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5. Probabilistic representation of Wiener Chaos solutions. As before, let 
F= (f2, {J~t}o<t<T,R) be a stochastic basis with the usual assumptions 
and let W = {wk(t),k > 1,0 < t < T} be a collection of standard Wiener 
processes on F. 

Consider the linear equation in R d 

(5.1) du = {aijDiDjU + b{DiU + cu + f)dt + {ai k DiU + v k u + g^) dw k 
under the following assumptions: 

(BO) All coefficients, free terms and the initial condition are nonrandom. 
(Bl) The functions aij = a>ij(t, x) are measurable and bounded in (t, x), and 
(i) 

\aij(t,x) - aij(t,y)\ < C\x - y\, x, y £ R d , < t < T, 

with C independent of t,x,y; 
(ii) the matrix (ay) is uniformly positive definite, that is, there exists 
a 5 > so that, for all vectors y G R d and all (t,x), a^yiyj > 5\y\ 2 . 
(B2) The functions bi = bi(t,x), c = c(t,x) and v k = v k (t,x) are measurable 

and bounded in (t,x). 
(B3) The functions ait = atk(t,x) are continuous and bounded in (t,x). 
(B4) There exists ap>rf+l so that the functions / = f(t,x) and g k = 
g k (t,x) belong to 

L p ((0,T);L 2 (R d )nL p (R d )). 

(B5) The initial condition u = u (x) belongs to L 2 (R d ) n W 2 (R d ), p>d+l, 
where W 2 is the Sobolev space {/ : /, Dtf, DiDjf G L p {R d )}. 

Under assumptions (B2)-(B4), there exists a sequence Q = {qk,k > 1} of 
positive numbers with the following properties: 

(PI) The matrix A with Aij = — (1/2) J2k>i Qk^ik^jk satisfies A^ (t, x)yiyj > 

0, x,y£R d ,0<t<T. 
(P2) There exists a number C > so that 

J2( su P\Qk^k(t,x)\ 2 + I ql\\g k \\ 2 L2(Rd) (t)dt\ <C. 

For the matrix A and each t, x, we have Aij(t, x) = ^^(t, x)aj k (t, x), where 
the functions are bounded. This representation might not be unique; 
see, for example, [9], Theorem III. 2. 2 or [36], Lemma 5.2.1. Given any such 
representation of A, consider the following backward Ito equation: 



X t , x ,i( s ) = x i + ( Bi(T,X t , x (T))dT + ( qk^ik(r,X ttX (T))d~Wk(T) 

Js k>l s 

+ a ik (T,X ttX (T))dtiT k (T); s£ (0,f),t€ (0, T], t-fixed, 
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where Bi = bi — J2k>i ^\ a ik v k an d w k ,k>l, are independent standard Wiener 
processes on F that are independent of w k ,k> 1. This equation might not 
have a strong solution, but does have weak, or martingale, solutions due to 
assumptions (B1)-(B3) and properties (PI) and (P2) of the sequence Q; this 
weak solution is unique in the sense of probability law ([36], Theorem 7.2.1). 
Let Q be the operator 

Q : {u a , a 6J}h {q a u a , a G J}. 

Theorem 5.1. Under assumptions (B0)-(B5), (5.1) has a unique Wiener 
Chaos solution u = u(t,x). IfQ is a sequence with properties (PI) and (P2), 
then u(i, •) G L2,q(W; L 2 (W i )), < t < T, and the following representation 
holds: 

u(t, x) = Q _1 E 1 f{s, X t , x (s))j(t, s, x) ds 

(5.3) +y~] [ q k gk(s,X t:X (s))j(t,s,x)d~uT k (s) 

k>i Jo 



+ u o (X tjX (0)) 1 (t,0,xm 

t < T, x G R rf , where Xf jX (s) is a weak solution of (5.2) and 
rt r t 



j(t,s,x) =exp( / c(T,X t>x (T))dT + ^2 f q k Vk(T,X tiX (T))a\vk(T) 
V s k>i Js 

(5-4) 

~lJ^<lkWk(r,Xt, x (r))\ 2 dT\ 



k>l 



Proof. It is enough to establish (5.3) when t = T. Consider the equation 
dU = (a ij D i D j U + biDiU + cU + f) dt 

( 5 - 5 ) +J2(<rikDiU + v k U + g k )q k dw k 

k>l 

with initial condition U(0,x) = uq(x). In Theorem 4.1, take Q = 1 and con- 
sider the normal triple (H%(R d ), L 2 (M. d ), H^iW 1 )). Then (5.5) has a unique 
Wiener Chaos solution U and 

U(t,-) €L 2 {W;L 2 (R d )), 0<t<T. 

By construction, process u = Q~ 1 U is the corresponding Wiener Chaos 
solution of (5.1). To establish representation (5.3), consider the function 
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Uh = E(U£(h)),h G 7i. According to Theorem 3.7, the function Uh is the 
unique solution of the equation 

dU h = (a^DiDjUh + biDiU h + cU h + f) dt 

(5.6) 

+ '52(<7ikD i U h + u k U h + g k )q k h k dt 
k>\ 

with initial condition Uh\t=o =uq. We also define 

Y(T,x) = [ T f(s,X T , x ( S )) 1 (T,s,x)ds 
Jo 

rT < 

(5.7) 9k(s,X Tx (s))j(T,s)q k dw k (s) 

k>i Jo 

+ u o (X T>x (0)) 1 (T,0,x). 

By direct computation, 

E(E(£(h)Y(T, x)\rf)) = E(£(h)Y(T, x)) = E'Y(T, x), 

where E' is the expectation with respect to the measure dF' T = £(h) cffiV and 
Pt is the restriction of P to J-jf. 

Under assumptions (B0)-(B5), the solution Uh of (5.6) is continuous in 
(t,x) and has a probabilistic representation via the Feynmann-Kac formula; 
see [21], Section 1.6. Using the Girsanov theorem ([14], Theorem 3.5.1), this 
representation can be written as Uh(T,x) = E'Y(T,x) or 

E(£(h)U(T, x)) = E(£(h)EY(t, x)\J^). 

By Remark 3.6, the last equality implies U(T,x) =E(Y(T,x)\J r ^ r ) as ele- 
ments of Zj2(W). 

Theorem 5.1 is proved. □ 

Example 5.2 (Krylov-Veretennikov formula). Consider the equation 

d 

(5.8) du = (ajjDiDjU + bjDju) dt + <jj k Diudw k , u(0,x) =Uq(x). 

k=l 

Assume (B0)-(B5) and suppose that a,ij(t,x) = ^<Ji k (t,x)<jj k (t,x). By Corol- 
lary 4.2, (5.8) has a square integrable solution such that 

E\\u\\l 2m (t)<C*\\u \\* 

By Theorem 4.1 with Q = 1, 

oo 

u(t,x) = ^2 ^2 U a (t,x)£ a 
n=l \ a \=n 
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u (x) 



(5.9) 




Pt,s n &jk n Dj 



• • • Ps 2 ,s 1 (Tik 1 DiP slfi u (x) dw kl (si) • • • dw kn {s n ) 



where P^ jS is the semigroup generated by the operator A = aijDiDjii+biDiU. 
On the other hand, in this case, Theorem 5.1 yields 



Thus, we have arrived at the Krylov-Veretennikov formula (cf. [20], Theo- 



6. Passive scalar in a Gaussian field. The following viscous transport 
equation is used to describe time evolution of a passive scalar 9 in a given 
velocity field v: 

(6.1) 9(t,x) = vA9(t,x) -v(t,x) ■ V0(t,x) + /(t,x); x£R d ,d>l. 

In the Kraichnan model of turbulent transport [10, 18], using the results 
from [1] and [23], (6.1) can be written as an Ito stochastic evolution equation: 

(6.2) d9(t, x) = (vA9(t, x) + f(t, x))dt- a k {x) ■ S79(t, x) dw k (t), 



(51) {wk(t),k > l,t > 0} is a collection of independent Wiener processes. 

(52) {a k ,k > 1} is an orthonormal basis in the space He, the reproducing 
kernel Hilbert space corresponding to the spatial covariance function C 



of v, so that a[(x)a{{y) = (»)?;%) = C l \x - y) and C^'(0) = 6%. 
The space H c is all or part of the Sobolev space H ( - d+ ^/ 2 (R d ;R d ), 
7 €(0,2). 




rem 4): 




•■•Ps 2 , Sl ^ DiU(p)(x) dw kl (si) ■ ■ • dw kn (s n ) . 



where 
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(S3) The initial condition 9q is nonrandom and belongs to L2(M d ). 

Using the Wiener Chaos approach, it is possible to consider velocity fields 
v that are even more turbulent, for example, 

(6.3) v i (t,x) = J24(x)Mt), 

k>0 

where {o~k,k > 1} is an orthonormal basis in L2(M rf ; M. d ). With v as in (6.3), 
the passive scalar equation (6.2) becomes 

(6.4) 9(t, x) = uA6(t, x) + f(t, x) - V6(t, x) ■ W{t, x), 

where W = W(t,x) is a (i-dimensional space-time white noise. 

In [6] and [34], a similar equation is studied using the white noise ap- 
proach. For more related results and references, see [12], Section 4.3. We 
will consider an even more general equation. 

Theorem 6.1. Suppose that v > is a real number, the functions c k {x) 
are bounded and measurable, 9q £ -^(R'O and f £ L2((0,T);ff 2 ~ 1 (M d )). 
Fix e > and let Q = {qk,k > 1} be a sequence so that, for all x, y £ ~R d , 

Mv? - ihli^ki^viyj ^ £ \y\ 2 - 

k>l 

Then there exists a unique Wiener Chaos solution of (6.2). This solution 
satisfies 

9eL 2 , Q (W; L 2 ( (0, T) ; E\ (R d ) ) ) n L 2 ,Q ( W; C ( (0, T) ; L 2 (R d ) ) ) . 

Proof. The result follows from Theorem 4.1, with A = uA, Mk = 
o-\Di, in the normal triple (JT^(R),L 2 (R),fl^ 1 (R)). □ 

Remark 6.2. If maxjsupj, |ct^.(x)| < Ck, k>l, then a possible choice 
of Q is q k = (Su) 1 / 2 /(d2 k C k ), < <5 < 2. 

If a k (x)o J k (x) <C a < +oo, i,j = 1, . . . , d, x £ M. d , then a possible choice of 
Q is q k = e{2v/{C G d)) 1 / 2 , < e < 1. 

When v = 0, (6.2) describes nonviscous transport [8] and can still be 
studied if interpreted in the Stratonovich sense: 

(6.5) du(t, x) = f(t, x) dt — o"fc(x) • V9(t, x) o dwk(t). 

In what follows, we assume that / = 0, each a k is divergence free and 
a k (x)a J k (x) =Sij. Then (6.5) has an equivalent ltd form 

(6.6) dB{t,x) = \Ae(t,x)dt-a k {x)D i 6{t,x)dw k {t). 

The following result (cf. [25]) summarizes the main facts about the Wiener 
Chaos solution of (6.6) 
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Theorem 6.3. In addition to (S1)-(S3), assume that each o k is diver- 
gence free. Then there exists a unique Wiener Chaos solution 6 = 6(t,x) of 

(6.6) . This solution has the following properties: 

(i) For eachtG [0,T], 0(t, •) E L 2 (R d ). 

(ii) For every (p £ C^(R d ) and all t £ [0,T], the equality 

(6.7) {6,<p)(t) = (0 ,<p) + ± f\e,Acp)(s)ds+ f\e,4D^)dw k (s) 

Jo Jo 

holds in L,2(W), where (•,•) is the inner product in L2(R d ). 

(iii) // 6 E W$(R d ), p>d+l, and X is a weak solution of 

(6.8) X tx = x+ [ a k (X S)X )dw k (s), 

Jo 

then 

(6.9) e(t,x)=E(9 (X t>x )\^), 0<t<T,xeR d . 

(iv) For 1 < p < oo and r£l, denote by L p ^(R d ) the space of measur- 
able functions with finite norm 




[ \f{x)\v{i + \x\ 2 y r / 2 dx. 

Ju d 



Then there exists a number K depending only on p and r so that 
( 6 - 10 ) nO\\i Mr)m (t) < e Kt \\e \\l p {r)[Rdy t > 0. 

In particular, if r = 0, then K = 0. 

Proof, (i) This follows from Theorem 4.1, because (6.6) is weakly 
parabolic and one can take Q = 1. 

(ii) Since for each k we have T)iO~\ = in the sense of distributions, the 
same arguments as in the proof of the second part of Theorem 3.4 result in 
(6.7). 

(iii) Equality (6.9) follows from Theorem 5.1 after observing that the 
time-homogeneity of a k allows us to rewrite the corresponding backward 
equation as (6.8). 

(iv) To prove (6.10), denote by St'.0o(-) t— >0(t, •), t > 0, the solution op- 
erator for (6.6). Direct calculations show that, for every r£f, assumptions 
(A1)-(A2) of Theorem 4.1 are satisfied with V = H 2 \ {r) (R d ) := {/ : /, Dtf £ 

L P)(r) (R d )}, H = L 2i(r) (M d ), V = H~} r) (R d ) and Q = l. Then S t is a bounded 

linear operator from L2^ r )(R d ) to L 2i ( r )(0 x an d 

\\ s t0o\\ L2t(r) {nxR d ) < e 2 *l|0o||L 2?Cr) (R<*) 
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for every r € M, where the number K 2 depends only on r, and 

ll/llL 2 ,(r)(^XK d ) :=E ll^llL,(r)(K d )' 

The calculations also show that K 2 = if r = 0. 

By (6.9), St is a bounded linear operator from 
and 

||'St^o||i 00 (nxR"*) - ll^o||L 00 (R d )- 
Interpreting the last inequality in the weighted spaces L^m with weight r = 
0, we interpolate between = £00,(0) and L 2 ^ r ) and derive (6.10) from the 
Stein interpolation theorem in weighted spaces (see, e.g., [2], Theorem 4.3.6). 
Theorem 6.3 is proved. □ 

Acknowledgment. We are indebted to the referee for a number of in- 
sightful suggestions. 
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